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Abstract 


A computational simulation of a transonic wind tunnel test section with 
longitudinally slotted walls is developed and described herein . The non- 
linear slot model includes dynamic pressure effects and a plenum pressure 
constraint , and each slot is treated individually . The solution is performed 
using a finite- difference method that solves an extended transonic small 
disturbance equation. The walls serve as the outer boundary conditions 
in the relaxation technique , and an interaction procedure is used at the 
slotted walls. Measured boundary pressures are not required to establish 
the wall conditions but are currently used to assess the accuracy of the 
simulation. This method can also calculate a free-air solution as well 
as solutions that employ the classical homogeneous wall conditions. The 
simulation is used to examine two commercial transport aircraft models 
at a supercritical Mach number for zero-lift and cruise conditions. Good 
agreement between measured and calculated wall pressures is obtained for 
the model geometries and flow conditions examined herein. Some local- 
ized disagreement is noted , which is attributed to improper simulation of 
viscous effects in the slots. 


Introduction 

The National Transonic Facility (NTF) at the 
Langley Research Center is designed to allow tran- 
sonic wind tunnel testing of model configurations at 
full-scale Reynolds numbers (ref. 1). Because full- 
scale Reynolds numbers can now be matched in test- 
ing, new data correction techniques have been devel- 
oped to further improve data quality. One correction 
addresses the interference that is induced by the test 
section walls. The ventilated test section has evolved 
in response to the requirements of transonic testing 
(ref. 2). A ventilated test section has walls that are 
either slotted or porous, and the entire test section is 
contained in a large, closed plenum chamber so that 
the tunnel can be pressurized. The NTF test section 
has six longitudinal slots in both the floor and ceil- 
ing that extend the full length of the test section and 
provide an open-area ratio of 6 percent of the wall 
area. The ventilated walls reduce but do not com- 
pletely eliminate wall interference in a test section. 
At transonic speeds, the wall interference problem is 
complicated by the possible nonlinear interaction of 
the model flow field with the near- wall flow field. 

Calculation of the flow over a slotted wall has tra- 
ditionally involved simplifications such as a homo- 
geneous wall treatment in which the flux through the 
open areas is distributed over the entire wall. Empir- 
ical correlations involving some features of the wall 
geometry such as openness ratio and slot width may 
be incorporated. A more complex treatment of slot- 
ted walls arises from a nonhomogeneous approach in 


which the flow through the slots is modeled in some 
fashion and the rest of the wall is treated as solid. 
With either approach, the boundary specification 
may rely on experimental data such as the pressure 
distribution on the walls. Because many accurate 
measurements are required for sufficient resolution of 
the three-dimensional flow T field, demanding experi- 
mental requirements are imposed. Wall interference 
approaches may be broadly classified as cither pre- 
dictive or assessment methods. Predictive methods 
are considered to be pretest examinations because 
measured data from the current test are not required 
for the simulation. Assessment methods may be 
viewed as post test examinations because current ex- 
perimental data are used in the expression of the wall 
boundary conditions. Several computational tech- 
niques have been used to examine wall interference ef- 
fects in the NTF (refs. 3-5). Reference 3 describes an 
assessment approach that uses measured wall pres- 
sure data and a homogeneous slot ted- wall assump- 
tion within a transonic small disturbance (TSD) flow 
solution. The method of reference 4 is a predictive 
technique that also combines a homogeneous outer 
boundary with a TSD solution. Reference 5 describes 
an assessment method in which a simulation of dis- 
crete slotted-walls is combined with a panel method 
formulation. 

Wall interference calculations may be performed 
by obtaining two computational solutions of the same 
aircraft configuration: one solution for an aircraft 
bounded by walls and the other for an aircraft in 


“free air.” Assuming that the pertinent flow physics 
are accurately modeled, the differences between the 
two solutions give an indication of the disturbances 
that the walls induce on the model flow field and 
provide insight for possible wall interference correc- 
tions to wind tunnel data. In the current work a 
predictive technique is developed in which a discrete 
slotted- wall model is incorporated within a TSD flow 
solution method. This approach features the ability 
to analyze aircraft flow fields with a nonlinear govern- 
ing equation that includes possible discrete-slot inter- 
actions and at the same time allows independence 
from experimentally measured wall data. A free-air 
outer boundary condition is included to enable tun- 
nel versus free- air comparisons to be made. The slot- 
ted walls may also be treated as homogeneous. The 
components of this method are described herein, and 
the effects of the parameters are demonstrated. Ulti- 
mately, the method is intended as a pretest indicator, 
showing whether wall interference effects are signif- 
icant for a given model under specified flow condi- 
tions. Present data are insufficient for determining 
whether the method can accurately predict interfer- 
ence effects on the aircraft model. In this work, the 
suitability of the wall boundary conditions for sim- 
ulation of the test section flow field is assessed by 
comparison of calculated results with available wall 
pressure data. 

Nomenclature 

a slot width (see fig. 3) 

b slot strip width (see fig. 3) 

Cp drag coefficient 

Ci lift coefficient 

Cm pitching-moment coefficient 

Cp pressure coefficient 

d distance between slot centerlines 

E unit step function (see discussion. of 

eq. (4)) 

h NTF test section half-width and half- 

height 

L homogeneous slot parameter (see 

eq. (3)) 

l slot depth (see fig- 3) 

M Mach number 

NTF National Transonic Facility 

n direction normal to wall 


Q normalized velocity potential of two- 

dimensional sink 

q slot volume flux 

R Reynolds number 

s transformed slot variable 

s p ,o value of s corresponding to y Pi0 

TSD transonic small disturbance 

U streamwise velocity 

V normalized velocity of two-dimensional 

sink 

x streamwise direction 

y spanwise direction 

yi vertical direction in slot local coordi- 

nate system 

y p plenum pressure surface location 

yp Q limiting value of yp 

z vertical direction 

Z[ spanwise direction in slot local coordi- 

nate system 

a angle of attack 

7 ratio of specific heats 

<5 normalized plenum pressure (see 

eq. (7)) 

£ vertical direction in computational 

domain 

T] spanwise direction in computational 

domain 

rja viscous reduction factor for slot width 

rj u viscous reduction factor for longitudi- 

nal velocity 

£ streamwise direction in computational 

domain 

p density ' 

a location of slot lip in transformed plane 

< p velocity potential 

(j) z ^ 0 spanwise velocity component at slot 

centerline 

(p velocity potential along slot centerline 

Subscripts: 

i inviscid value at slot entrance 

l local slot coordinate system 


2 


n 

o 

P 

S 

w 

w 

x 

y 

z 

oo 


derivative with respect to normal 
direction 

value at slot origin 
value in plenum 

representative slot value, including 
viscous effects 

value at w r all 

value at wing 

derivative with respect to x direction 
derivative with respect to y direction 
derivative with respect to z direction 
value at free-air or reference condition 


Procedure 

The code developed herein consists of a tran- 
sonic small disturbance (TSD) aerodynamic analysis 
method that allows different outer boundary treat- 
ments to simulate both free-air and test section flow 
over an aircraft configuration. A discrete-slot flow 
model is incorporated to allow prediction of the be- 
havior of the flow near the slotted walls of the NTF 
test section (fig. 1). In a test section simulation, the 
TSD and slot flow solutions are performed simultane- 
ously using an interaction technique. The TSD code 
is used to solve the aircraft flow field and supply the 
slot flow calculation with normal volume flux values 
at the wall. The slot geometry and normal flux are 
used by the slot flow model to determine the potential 
along each slot. These values of potential are then 
used by the analysis code as an updated outer bound- 
ary condition along the slots. On the solid portions 
of the top and bottom walls and on the side wall, the 
standard tangential-flow boundary condition is im- 
posed. Thus, in a test section simulation, the outer 
boundary condition for the aerodynamic analysis is of 
mixed (nonhomogeneous) character, with potential 
specified in the slot regions and the normal gradient 
of potential specified on the solid portions. 

In this section, descriptions of the different ele- 
ments of the wall interference prediction code are 
given. The section briefly reviews the classical homo- 
geneous wall boundary conditions included in the 
code, then proceeds with a development of the 
discrete-slot flow model. Next, the focus shifts to 
the aerodynamic analysis code that is used. A brief 
outline of the formulation is given, followed by the 
modifications and additions necessary for simulating 
the test section flow. Finally, the interaction pro- 
cedure between the wall solution and the flow field 
solution is discussed. 


Wall Boundary Conditions 

Homogeneous wall boundary conditions . 
Several forms of the classical linear homogeneous 
wall boundary conditions are commonly used. In the 
present work, as in reference 4, the integrated, gen- 
eralized linear homogeneous formulations given by 
Keller in reference 6 are used. At a solid bound- 
ary, the condition to be satisfied is that no normal 
flow passes through the boundary. Thus, the nor- 
mal velocity component vanishes at the wall, which 
is mathematically stated as 



(i) 


At an open-jet boundary, no pressure gradient is al- 
lowed to exist across the boundary. In the linear ap- 
proximation, the pressure coefficient is proportional 
to the streamwise gradient of the velocity potential 
( C p ss —2 4> x ). Because the external flow is quies- 
cent ( C p = 0), the approximation to this condition 
is given by 


At a wall with several longitudinal slots, the lin- 
earized homogeneous condition states that the pres- 
sure gradient is balanced by the streamline curvature. 
The expression of this condition is 


4>x T L 


< 9 2 (/> 
dn dx 


= 0 


where L is a geometry-dependent slot 
given by 



parameter 


The integrated form of this condition as applied in 
the present work is written as 


deb 

on 


( 3 ) 


Discrete-slot boundary condition . The 

discrete-slot model was developed over a period of 
years by several researchers, and the earliest pre- 
sentation was a 1975 paper by Berndt and Sorensen 
(ref. 7). The effort evolved from what was initially 
an attempt to develop an improved homogeneous 
wall condition using a theoretical approach that 
wus consistent with certain observations of experi- 
mental slot flow. Berndt developed an inviscid, three- 
dimensional theoretical representation of the flow 
through longitudinal slots in which the flow leaving 
the test section through each slot is treated as a thin 
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jet while the return flow admits quiescent plenum air 
as shown in figure 2 (ref. 8). The model incorporates 
a quadratic crossflow velocity term (normal to the 
slotted wall) that is neglected in the linearized theo- 
ries. Karlsson and Sedin (ref. 9) presented an imple- 
mentation of the Berndt theory for an axisymmetric 
test section in which the steps in the solution process 
were detailed, although the research was directed to- 
ward the inverse problem of designing interference- 
free slot shapes. A formulation that includes a 
simple accounting for viscous effects in the axi- 
symmetric model is presented in reference 10. In a 
separate effort , the basic ideas of the Berndt inviscid 
slot flow model are extended to a formulation that 
does not require an axisymmetric test section geo- 
metry (ref. 11). For the axisymmetric case, Berndt 
performed a matched asymptotic expansion of the 
test section and slot flow fields that resulted in an im- 
proved homogeneous outer boundary condition. In 
the extension to nonaxisymmetric test sections, a 
homogeneous outer condition no longer results. The 
asymptotic matching is performed in a more local- 
ized sense; each slot is treated separately and has its 
own local matching. The effect is that each slot is 
“smeared” somewhat in the lateral direction. In dis- 
tinct slot regions that are wider than the slots, flow 
normal to the wall is permitted. The slot regions 
are separated by solid regions in which there is no 
flow normal to the wall. The result is a mixed outer 
boundary condition for the test section flow calcu- 
lation in which is specified on the slot strips and 
<f> n is specified in the solid- wall regions (fig. 3). The 
viscous parameters of reference 10 are incorporated 
in the nonaxisymmetric slotted-wall model in refer- 
ence 12. The slot model of references 11 and 12 is 
used in the current work. 

In formulating the slot boundary condition, the 
intent is to create a velocity profile across the slot 
strip that simulates the flow through a slot in a far- 
field sense (as seen from the model location in the 
test section). This approximation is carried out by 
defining a slot strip that is slightly wider than the 
actual slot (fig. 3). A solution for the flow through 
a slot involves a solution of two first-order ordinary 
differential equations along the slot centerplane. One 
equation gives the location of a surface along which 
the plenum pressure applies. This surface, denoted 
y p , may be considered as the boundary between the 
test section flow and the plenum (fig. 2). The plenum 
is assumed to be quiescent; no modeling of the flow is 
done inside the plenum, and the value of the plenum 
pressure is an input constant that is imposed on the 
plenum side of y p . The limiting depth that y p is 
allowed to extend into the plenum is denoted by 


a similar surface y Py0 , which is a function of the 
slot geometry only. When y p reaches y }KO , the slot 
flow is assumed to separate into the plenum as a 
thin jet (fig. 2). The second differential equation 
is a normal momentum equation that expresses the 
pressure difference maintained by the slot. The 
momentum equation is applied through the slot from 
the outer region of the test section to the plenum 
pressure surface. With a known plenum pressure and 
the momentum equation, which gives the pressure 
gradient maintained by the slot, the pressure at 
the boundary of the test section (that is, at the 
slot entrance) is determined. Because the pressure 
coefficient is cast in terms of the gradient of potential, 
the integrated solution actually yields the potential 
along the slot centerline y>. After this centerline 
potential is known, the potential across the entire 
slot strip is given by the asymptotic matching. 

The details of the slot model are now presented. 
Each slot is treated separately and has its own local 
coordinate system in which x corresponds to the free- 
stream flow direction, yi to the normal or crossflow 
direction (positive outward into the plenum), and zi 
to the spanwise direction (fig. 3). For a given slot, the 
wall surface is located at y/ = 0, the slot centerplane 
is given by Z[ = 0, and Z[ ranges between +6/2 and 
—6/2, where b is the width of the slot strip. The 
slot flow formulation in the slot centerplane is based 
on a conformal transformation to a potential flow 
source/sink in the far-field limit (infinitely deep in 
the slot). (See ref. 8.) Let Q denote the normalized 
velocity potential of an isolated source /sink of unit 
width and flux, and V denote the corresponding 
velocity. The equation that describes the location 
of the plenum pressure surface is (ref. 11) 

dy P q s V ^ V P da s 

= E — {V 

ax a s a s ax 

where a is the slot width and q is the normal volume 
flux per unit length through the slot. The subscript 
s denotes a representative slot value in which viscous 
effects have been included. The viscous model will be 
discussed below. Equation (4) states that the slope 
of the plenum pressure surface is proportional to the 
local flow angle. The second term on the right side is 
a one- dimensional correction for varying slot width, 
where E = 1 inside the slot and E = 0 outside the 
slot. The expression defining y PyQ (ref. 9) is 


Vyp 

a s 


l , ln(2) 
a s 7r 


1 + 2 exp 


7 tI 

1 + — 
0 * 


( 5 ) 


where l is the depth of the slot. When the value of y p 
given by equation (4) exceeds y P)0 , then y p is replaced 
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by y p ^ 0 . Thus, an integration of equation (4) along 
the slot, subject to the constraint of equation (5), 
yields the location of the plenum pressure surface. 

The crossflow momentum equation in the slot 
centerplane is developed by considering a two- 
dimensional crossflow plane at a given x location. 
The equation (ref. 11) is 


d(v-(p 0 ) p s U s d f U s f 
dx p a 0 Uoc dx U. x 



5 — f^) 2 f — - ^ — ^) 2 -« 

* Poc V Uog / V Q-s a s dx J 


( 6 ) 


where ip Q is the value of ip at the slot origin. The 
normalized plenum pressure parameter 6 is defined 
as 


Pp- Poo 
Poc^oc 


( 7 ) 


and is equal to one half of the pressure coefficient 
as typically defined. Also, note the presence of the 
quadratic crossflow velocity term on the right side of 
equation (6). Equation (6) is integrated along the 
slot to give the distribution of the slot centerline 
potential <p. The asymptotic matching defines the 
relationship between ip and the potential across the 
entire slot strip. This expression is written in the 
local slot coordinate system as (ref. 11) 


where U is the local longitudinal velocity component 
at the slot entrance, a is the slot width, and the sub- 
script i denotes the inviscid value at the slot entrance. 
To keep the slot model as simple as possible, the vis- 
cous reduction factors are assumed to have constant 
values along the slots. 

Small perturbation approximations are used to 
derive expressions for the necessary velocity and 
density ratios along the centerline of each slot. Using 
small perturbation approximations, the local density 
and streamwise velocity are 

— = 1 (11) 

Poo 

Vi 

+ (12) 

'-'OC 

Using equation (9) and assuming constant total 
enthalpy and pressure, the characteristic slot flow 
density can be written as (ref. 10) 


Ps 

Pi 




(13) 


where 

M{ = A/qc [1 + (7 — 1) 0 X ] 1/2 (14) 

Combination of equations (11) and (13) yields the re- 
quired expression for the slot density ratio, and the 
streamwise component of the local slot flow veloc- 
ity is written using equations (9) and (12). These 
expressions are 


J> M z l ) = -t?±( Z j)\ ip + ( fi ZtOZl ( 8 ) 

where (j> z ^ 0 is the sidewash (spanwise) velocity. The 
sidewash velocity is calculated along the centerline of 
each slot using the potential from the TSD solution. 
The values of potential given by equation (8) are ul- 
timately imposed along the slot strips as a boundary 
condition for the TSD solution. 

Several physical phenomena are considered in ex- 
tending the slot flow model to include viscous effects 
(ref. 10). The streamwise velocity at the slot entrance 
is assumed to behave as in a shear flow, with a corre- 
sponding loss of longitudinal momentum. Similarly, 
the crossflow velocity through the slot is assumed to 
have a profile similar to that of a boundary layer flow, 
with an associated increase in core velocity. These 
behaviors are simulated by viscous reduction factors 
for the longitudinal velocity 1 } u and slot width rj a , 
defined as 

U s = T) u Ui (9) 

( 10 ) 


Ps l - M^4> x 

(15) 


u s 


Yj (1 4" ^x) Vu 

u oo 

(16) 


The normal volume flux q is required in the 
solution of the slot flow equations. Estimates for q 
are provided through spanwise integration, across the 
slot strip, of the normal velocity component at the 
wall, or 

Pi Uoo f + b / 2 

qs = n 7T / ^ dZ l ( 17 ) 

Ps Us J-b/2 

The normal velocities over each slot strip are ob- 
tained from the TSD solution of the test section flow 
field. In reference 11 researchers found that using nu- 
merically integrated values for q in this manner gave 
results that were insensitive to the value chosen for b. 

The slot flow formulation is completed with the 
addition of the relations describing the source/sink 
potential and velocity functions Q and V. As men- 
tioned above, these functions are derived through use 


5 


of a conformal transformation. An outline of the for- 
mulation is given in the appendix of reference 8, and 
the results are summarized here for convenience. In 
terms of the transform variable 5 , the source/sink 
functions (ref. 8) are given by 

Q — — a In s (18) 


V = 


a 

\/s 2 + G 1 


(19) 


The parameter <j denotes the location of the corner of 
the slot lip in the transformed plane, where a — l/n. 
Finally, the relationship between the physical and 
transform variables (ref. 9) is expressed by 


yp = 

0'S 


\/s 2 + <t 2 + a In 



(20) 


To summarize, a simultaneous solution of equa- 
tions (4) and (6) yields the distribution of <p along 
the centerline of a slot. Equations (18) and (19) and 
the auxiliary expression given by equation (20) define 
the functions that appear in these differential equa- 
tions. After <p is determined, the potential over the 
entire slot strip is specified by equation (8). A discus- 
sion of the integration technique used in the current 
work is presented in appendix A. 


Aerodynamic Analysis 

The aerodynamic analysis calculation is per- 
formed with a three-dimensional TSD method devel- 
oped by Rosen (ref. 13) as an extension of the work 
of Boppe (ref. 14). The method is based on a formu- 
lation of the TSD equation that has been extended 
by the retention of some higher order terms usually 
neglected in TSD formulations. This method allows 
the treatment of complex configurations that include 
a wing, fuselage, pylons, and underwing stores. A 
global coarse grid extends over the entire flow do- 
main, and much finer embedded grids are used to 
give good resolution in the wing, pylon, and store re- 
gions. In the coordinate system used, x denotes the 
direction of the free stream, y is the spanwise direc- 
tion, and z is the vertical direction. The governing 
equation is written as 

T(j> XX + U(f>xy + V (j)yy + <f>ZZ = 0 ( 21 ) 


where 

T = 1 - Ml, - (7 + 1) Ml ,4> x 
U = -2Ml4>y 


V = 1 - ( 7 - 1 ) M^<p x 

The terms containing < px<t>xx* fiyfixy, and fix^yy are 
the higher order terms that are typically neglected 
in TSD formulations. The term improves the 

transition between subsonic and supersonic zones, 
and the (pyCpxy and 4>x4>yy terms are retained to aid 
in the resolution of swept shock waves. An upwind, 
rotated-differencing scheme is used in the numerical 
relaxation procedure to ensure that the proper do- 
main of dependence is maintained in supersonic flow 
regions. Rosen developed a scheme specifically for 
this modified TSD formulation in which the rotation 
is determined using the coefficients T, [/, and V. 
As a result of the * extended TSD formulation and 
the rotated-differencing scheme, the method can ac- 
curately treat highly swept and tapered wings at 
high transonic Mach numbers with good numerical 
stability and accuracy (ref. 13). 

Wall Interference Code 

The code created in the course of the present work 
is constructed by making additions and modifications 
to the Rosen code. These alterations include a 
new grid generation scheme for the global coarse 
grid, the calculation and imposition of the different 
outer boundary conditions, the simulation of the 
model support sting, and the ability to simulate a 
finite-length test section. Each of these features is 
discussed here. 

Grid generation . Several factors must be con- 
sidered in developing the grid generation scheme for 
the global coarse grid. To examine wall interfer- 
ence effects, solutions obtained from the test section 
simulation are compared with free-air solutions. To 
eliminate grid-biasing differences from these compar- 
isons, the test section grid is wholly contained within 
a larger free-air grid with outer boundaries that tend 
toward infinity (ref. 4). Grid planes are located at 
each boundary of the test section, including the test 
section walls as well as appropriate inflow and out- 
flow locations. At the boundary of the test section 
grid, a smooth transition to the free-air grid is re- 
quired. The total number of grid points should be 
kept to a reasonable number, while maintaining ade- 
quate grid resolution at the walls without sacrificing 
grid resolution at the model location. The spacing 
in the original Rosen code was used as a guide to 
determine the grid spacing near the model in each 
of the coordinate directions. In both the Rosen code 
and the present code, a vertical symmetry plane is 
assumed to exist at the centerline of the aircraft. As 
a result, only half the flow domain is modeled. A 
brief discussion of the grid scheme will now be given, 
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and details of the grid generation are contained in 
appendix B. 

The grid in the streamwise (x) direction is divided 
into five regions and is symmetric about the midpoint 
of the central region. (See fig. Bl.) The central 
region, which contains the wing, consists of uniform 
spacing that results in the finest resolution of the 
streamwise grid. The two regions that bound the 
central region extend from the respective edges of the 
central region to the test section inflow and outflow 
planes. A mild stretching function is used in these 
regions so that the spacing is still reasonably fine 
near the test section boundaries. In the regions 
outside the test section, a strong stretching function 
is used to rapidly extend the grid toward upstream 
and downstream infinity. For the aircraft geometries 
currently being examined, a total of 111 grid points 
are used in the streamwise direction for a free- air 
calculation. Of these grid points, 17 are clustered 
within the wing region. For a test section calculation, 
a total of 84 streamwise grid points are used with the 
current geometries. 

A discussion of the spanwise grid generation re- 
quires a brief explanation. In the original Rosen 
scheme, the spanwise spacing is determined by the 
wing span, with the geometric wingtip midway be- 
tween two grid planes. No attempt is made to pre- 
dict the actual flow at the tip since the region is dom- 
inated by rotational and viscous effects that cannot 
be modeled with a TSD approach. In the context 
of potential flow, the intent is to avoid performing a 
calculation too near the tip while presuming that the 
lift goes to 0 at the tip. Thus, location of the wingtip 
midway between two grid planes gives the best treat- 
ment of the tip region. In the current scheme con- 
venience dictates that the spanwise grid spacing is 
determined by the test section geometry rather than 
by the wing geometry. The geometry of the NTF 
test section is such that slots are located at one-sixth, 
one-half, and five-sixths of the test section half-width 
(fig. 1). Thus, a uniform grid spacing with an incre- 
ment equal to some fraction of one sixth of the test 
section half-width ensures that a grid plane is aligned 
with each slot as well as with the test section side- 
wall. A consequence of the present technique is that 
the wingtip has an arbitrary location between two 
grid planes. Care must be taken to treat the tip re- 
gion adequately, as discussed above. To generate the 
spanwise grid, the spacing in the test section region 
is varied by altering the total number of grid planes 
in increments of six until enough planes fall on the 
wing (typically 15-20). To avoid performing an in- 
viscid calculation too near the wingtip, the outer- 
most of these wing grid planes may be ignored if too 


close to the geometric wingtip (that is, within 2 per- 
cent of the semispan). The lift is assumed to go to 0 
at the tip, as in the original method. This proce- 
dure was found to give spanwise lift results similar 
to those of the original method. Outside the wall, a 
strong stretching function is used again to extend the 
grid toward infinity. For the free-air calculation, a to- 
tal of 45 grid points is used in the spanwise direction. 
Of these points, 20 fall on the model (18 on the wing) 
for the geometries currently examined. When a test 
section simulation is performed, the total number of 
spanwise grid points decreases to 37. 

Spacing grid points uniformly in the vertical di- 
rection while maintaining adequate resolution near 
the wing would require an excessive number of grid 
points (that is, double the number used in the span- 
wise direction for the same spacing). Ideally, the grid 
would have the finest spacing near the model, fine 
spacing near the wall, and coarser spacing between 
these regions. In the grid generation scheme that was 
developed, three zones are used inside the test sec- 
tion to give the desired spacing (fig. Bl). The central 
zone is symmetric about the wing plane and extends 
over half the test section height. A mild stretching is 
used in this region, which gives the finest spacing at 
the wing. Between the edges of the central zone and 
the walls, a similar stretching is used to slightly de- 
crease the spacing as the walls are approached. The 
stretching functions in the two regions (near the up- 
per and lower walls) may differ slightly because the 
wing plane may not be in the exact center of the test 
section. Outside the walls a strong stretching func- 
tion is used to give the required extension toward in- 
finity. The free-air grid currently used has 63 points 
in the vertical direction. With the test section grid 
the number of points decreases to 51. 

The above discussion refers to the work that 
was done on the global coarse grid. Note that an 
embedded fine grid is used to resolve the details of 
the wing flow field (ref. 13). The wing fine grid is 
automatically generated, based on the wing geometry 
and the global coarse grid. Fine grid planes coincide 
with the spanwise coarse grid planes that lie on 
the wing; thus, for the current geometries there are 
18 spanwise fine grid planes. Each fine grid plane 
has 130 points in the streamwise direction. Of these, 
100 of the points are located along the chord at 
equally spaced increments. The extent of the fine 
grid in the streamwise and vertical directions is a 
set percentage of the local wing chord length. The 
fine grid extends from 20 percent chord upstream of 
the leading edge to 10 percent chord downstream of 
the trailing edge. In the vertical direction, the fine 
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grid extends from 10 percent chord below the wing 
to 30 percent chord above the wing. 

Test section features . Two features added to 
the code to approximate the flow in a test section 
more accurately arc simulations of a model support 
sting and a finite-length test section. The cross- 
sectional area of the sting is generally of the same 
order of magnitude as that of the model fuselage. 
Thus, in the test section environment the sting in- 
duces a significant blockage downstream of the model 
and influences the wall pressure coefficients. A simu- 
lation of the sting is included here because the results 
of the code are compared with measured wall pres- 
sures. In the code, the sting is treated like an axi- 
symmetric body that is an extension of the fuselage, 
which is an adequate model because the fuselage it- 
self is treated as an axisymmctric body. Note that 
this representation is realistic for configurations such 
as the Pathfinder I model with its rear-mounted axi- 
symmetric sting (fig. 4(a)). With blade sting ar- 
rangements such as the Boeing 767 model and its 
high- mounted sting (fig. 4(b)), the flow in the sting 
region will not be accurately simulated. At a min- 
imum the present simple sting representation is 
expected to approximate the gross blockage effects. 

In the test section simulation, the test section in- 
flow plane coincides with the beginning of the slots 
and the outflow plane with the end of the slots. From 
a numerical standpoint, these choices are found to 
give the best behavior in the inflow and outflow re- 
gions. To justify these choices in a physical context, 
the geometry of the NTF test section must be con- 
sidered. In the NTF test section, the slots begin 
6 in. downstream of the test section entrance and 
gradually open from an initial width of 0 (table 1). 
Because the slots begin so close to the test section 
entrance, the slot beginning is a suitable location for 
the computational inflow plane. The slots terminate 
in a reentry region in which a step discontinuity in 
wall height is coupled with reentry flaps to guide in- 
flow from the plenum (fig. 1). The reentry region is 
dominated by viscous mixing near the wall, especially 
downstream of the step. Because accurate modeling 
of the viscous mixing region is difficult in an inviscid 
calculation, this region is excluded from the calcu- 
lation for convenience. The calculated model forces 
and moments were insensitive to the movement of the 
computational outflow plane location over the extent 
of the reentry region, and the calculated wall pres- 
sure coefficients showed only a localized sensitivity to 
the outflow location. Thus, the outflow plane may be 
located at the end of the slots, just upstream of the 
test section reentry region. A corresponding outflow 
boundary condition is discussed below. 


Imposition of boundary conditions . For a 
free-air calculation, the outer boundary conditions 
are the same as those in the original code (ref. 13). 
At all outer boundaries other than downstream in- 
finity, the boundary condition is no perturbation po- 
tential, or (j> — 0. At downstream infinity, the condi- 
tion is <f>x = 4>xx — <!>xy = 0. The imposition of this 
condition means that a Laplace equation is solved in 
the outflow^ plane (that is, (j) yy -j- <p zz = 0) to allow 
for the wing wake and circulation. At the centerline 
symmetry plane (y = 0), the condition is that the 
spanwise flow component vanishes, or (j) y = <p xy = 0. 
The test section simulation requires modifications in 
the formulation and imposition of the outer bound- 
ary conditions. At the inflow plane, the condition 
is still 0 = 0, and the symmetry plane condition is 
also unchanged. A wall boundary condition is im- 
posed at each of the walls. The capability exists in 
the code to treat each wall using the homogeneous 
solid, open, or slotted boundary conditions and ad- 
ditionally, to treat the floor and ceiling using the 
discrete-slot model. The homogeneous conditions are 
applied in the same manner as in reference 4. When 
the discrete-slot model is used, the outer boundary 
condition becomes mixed in character. In the solid- 
walled portions, the homogeneous solid-wall condi- 
tion (<fi n = 0) is applied. In the slot regions, the 
value of potential obtained from the slot solution is 
imposed. 

The treatment of the out flow boundary in the con- 
text of a finite-length test section is more difficult. In 
terminating the solution ahead of the reentry region, 
the flow is not uniform at the outflow boundary. Ex- 
amination of wall pressure data from the NTF shows 
that the wall pressure coefficient in the region of the 
outflow plane is still changing; that is, in general <j> xx 
is nonzero. In the current work, the condition im- 
posed is that 0 XX is equal to some constant in the 
outflow plane. The value of the constant is deter- 
mined by comparison with experimental wall pres- 
sure data. The value of the constant affects only the 
flow near the outflow plane and not the model forces 
and moments. Note that 0 XX represents a longitu- 
dinal flow acceleration, approximating the derivative 
of the pressure coefficient with respect to z. Con- 
ceivably, the proper specification of this slope in the 
outflow plane is a realistic alternative to modeling the 
viscous reentry region because it mimics the overall 
effect of the reentry region on the flow. Although 
<f) xx is unlikely to be constant over the entire outflow 
plane, such a simple condition may be sufficient to 
characterize the local flow behavior. 

Consideration of plenum pressure. The 
value of the plenum pressure must be specified in 
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the slot flow model. The plenum pressure governs 
the relative amounts of outflow and inflow through 
the slots as well as the total amount of flux through 
the slots. In reality the plenum pressure evolves pas- 
sively to a value that corresponds to an equilibrium 
in the slot and reentry flow. The plenum pressure 
depends on several factors such as the free-stream 
test conditions, the model geometry, the settings of 
the test section walls and reentry flaps, and the vis- 
cous mixing in the plenum and the reentry region. 
Thus, theoretical calculation of the plenum pressure 
would be difficult. In the inviscid method of refer- 
ence 5, a unique value of plenum pressure resulted in 
the proper numerical behavior of the slot fluxes for 
a given flow configuration. Because the solution was 
sensitive to the value of the plenum pressure, an in- 
verse approach was used in reference 5 in which the 
slot flux and the reentry-region streamline slope were 
specified, and the plenum pressure was determined 
in the course of solution. In the present method, the 
intent is to simulate the effect of the reentry region 
by specifying the outflow value of <f> X x , as discussed 
above, and to allow the slot fluxes to evolve, based on 
the specified plenum pressure. Because the plenum 
pressure is used to artificially force the proper slot 
flow behavior, including the effects of viscosity, the 
specified value of the plenum C p is generally not the 
same as the experimental value. 

The specified plenum pressure is felt through- 
out the test section. This effect is illustrated by 
consideration of the linearized expression for the 
pressure coefficient, C p = — In the linearized 
sense a nonzero plenum C p induces a longitudinal 
velocity perturbation <j> x through the test section. 
This perturbation results in an increment to the 
Mach number such that M is no longer equal to 
A/oc at the inflow plane. The free-stream Mach 
number must be adjusted to ensure that the Mach 
number at the test section entrance (that is, the 
test Mach number M ) remains the same when the 
plenum C p is changed. These Mach number adjust- 
ments can be determined using an expression derived 
from the isentropic pressure relation (ref. 15). The 
appropriate equation is 

A/ 2 = [(l + ^Ml) ( 7 A O + 1) (1_7>/7 - l] ^ 

( 22 ) 

Interaction and solution . The solution pro- 
cedure for a test section simulation requires an 
interaction technique: the aerodynamic solution and 
the slot flow solutions are coupled to obtain an over- 
all solution. The aerodynamic solution is an itera- 
tive procedure, while the slot flow solution involves 


a marching procedure from known conditions at the 
slot beginning. The slot flow and the test section flow 
are coupled through the wall values of the potential 
and the normal gradient of potential. A slot flow so- 
lution requires specification of the normal gradient 
of potential along the slot. The slot solution then 
provides the value of potential along the slot, which 
is imposed as an updated outer boundary condition 
on the aerodynamic solution. This iterative inter- 
action procedure continues with an exchange of up- 
dated values of potential and the gradient of potential 
until the aerodynamic solution converges. The wall 
pressure coefficients are calculated as the final step of 
the solution. (With a free-air calculation, the pres- 
sure coefficients arc calculated at locations equivalent 
to the wall positions.) The typical expression for the 
small perturbation pressure coefficient is used. The 
appropriate expression is 

Cp = -[2 cj> x + (1 - a4) 4 + 4 + 4 } (23) 

An interaction step begins with a calculation of 
the normal potential gradient (j> n over each slot- 
ted wall. At each longitudinal (x) grid station 
of the global grid, a span wise cubic spline inte- 
gration of the (j) n values is performed across each 
slot strip to obtain the local values of the slot vol- 
ume flux q. This curve-fit integration is used be- 
cause the edges of the slot strips do not neces- 
sarily coincide with spanwise grid-plane locations. 
The result is the distribution of q along the center- 
line of each slot with q values at the global coarse 
grid x locations. The slot flow equations are 
integrated using a fourth-order Runge-Kutta pro- 
cedure that requires the functional values to be given 
at uniform increments. Thus, a cubic spline interpo- 
lation of the q values is performed along each slot 
to give values of q at the proper intervals. The slot 
equations are then integrated with the Runge-Kutta 
method to yield values of <p along each slot. A cubic 
spline interpolation is performed to yield values of <p 
at the global coarse grid x locations. Equation (8) is 
then applied with these values to give the potential 
over each slot strip. The imposition of the potential 
values along the slot strips as a boundary condition in 
the TSD code requires a numerical underrelaxation 
to achieve a stable procedure. The value of the relax- 
ation factor depends on the frequency with which the 
slot solution is updated. An updated slot solution is 
currently performed with every other iteration of the 
TSD code to reduce the time required for the slot cal- 
culation by half. Any further reduction in the update 
frequency adversely affects the solution convergence. 
With an update frequency of every second iteration, 
a slot strip relaxation factor of 0.20 is currently used. 
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Because the slots begin at the inflow plane, the 
initial conditions for each slot integration are that 
the plenum pressure surface is level with the wall 
{yp — 0) and that the potential has a value of 0 
((f) = (p 0 = 0). Care must be taken in beginning 
the slot integration. The NTF slot geometry is 
such that the slots gradually open from an initial 
width of 0. Because the slot width is used as a 
normalization factor in the slot flow equations, a 
singularity exists at the slot origin. In the present 
method, the slot solution actually begins slightly 
downstream of the geometric slot origin so that the 
solution is initiated with a nonzero slot width. This 
approach is analogous to allowing no flux through 
the slot when the slot width is 0. The step size 
chosen for the Runge-Kutta integration is such that 
228 intervals are used along the slot. This number is 
judged to be sufficient because doubling the number 
of intervals results in no changes in the solution. 

The code has core memory requirements of ap- 
proximately 2.1 megawords and is presently im- 
plemented on a Cray Y-MP computer at Langley 
Research Center. The code is allowed to run until the 
prescribed number of iterations arc completed. For 
the results shown here 100 crude and 700 crude/fine 
iterations are used. This number of iterations is 
sufficient to allow the maximum change in poten- 
tial to decrease by at least three orders of magni- 
tude. The solution time varies, depending on the 
type of solution, with typical solution times on the 
Cray Y-MP as follows. An inviscid solution with dis- 
crete slotted walls requires 490 CPU sec. An inviscid 
calculation with homogeneous wall boundary condi- 
tions takes 440 CPU sec, which indicates that the 
discrete-slot calculations increase the solution time 
by about 11 percent. An inviscid free-air calculation 
requires 690 CPU sec. This 57-percent increase over 
the other homogeneous conditions is due to the addi- 
tional number of grid points required by the free-air 
solution. 

Analysis 

An analysis of the test section simulation is con- 
ducted in this section. Elements of both the test 
section and the aircraft model representations are 
considered, with the goal of comparing the simula- 
tion results with measured data. Test section wall 
pressure coefficients arc the primary data used for 
comparison. Note that a limited amount of suitable 
NTF wall pressure data currently exists because only 
simple transport aircraft configurations (fuselage and 
wing) are considered for calculation in this study and 
because acceptable tare pressure data must also ex- 
ist (as discussed below). The data from single tests 


of the Pathfinder I and Boeing 767 models are used 
here (fig. 4). To compare measured and calculated 
wall pressures, the angle of attack in the simulation 
is adjusted to obtain a calculated Cl that matches 
the measured Cl . Early work verified that the calcu- 
lated wall pressures are largely a function of the lift 
generated by the model. Although the code includes 
the provision to calculate the wing boundary layer, 
the boundary layer affects the wall C v values only 
indirectly through a slight change in the wing lift. 
Because the lift is easily adjusted by changing the an- 
gle of attack, the code is executed with no boundary 
layer calculation for the purpose of these studies. Ad- 
ditionally, note that the model configuration in the 
Pathfinder I test includes a horizontal stabilizer. Be- 
cause this method is capable of treating only a single 
lifting surface, the effect of the horizontal stabilizer 
cannot be modeled. 

This section begins with a brief comparison of the 
results of using different homogeneous wall boundary 
conditions, then presents a discussion of the mea- 
sured wall pressures that arc used herein for vali- 
dation. The behavior of the test section simulation 
is then examined in detail with a demonstration of 
the effects of varying simulation parameters on the 
solution. Among these parameters are components 
of the discrete-slot flow model and elements of the 
test section representation. Next, the calculated and 
experimental pressure coefficients are used to select 
the combination of model parameters that gives the 
best simulation of the test section flow field. Finally, 
test section and free-air calculations are compared 
to illustrate the influence of the test section on the 
aircraft flow field. 

Homogeneous Boundary Calculations 

The results of using various homogeneous wall 
boundary conditions are shown to illustrate the gen- 
eral influence of the outer boundary on the solution 
and to check in a simple way the behavior of the simu- 
lation method. The calculations arc performed with 
the Pathfinder I model at a Mach number of 0.80 
and an angle of attack of 2.70° to approximate cruise 
conditions for this aircraft. In figure 5, the wall pres- 
sure coefficients with the open-jet, solid, and slotted 
boundary conditions are compared with those calcu- 
lated at locations equivalent to the walls using the 
free- air boundary condition. The sidewall is treated 
as solid in each of the test section calculations. The 
fuselage of the Pathfinder I model is approximately 
between stations 11 and 15, and the wing is cen- 
tered near station 13. The figure shows that solid 
boundaries result in the greatest flow acceleration 
through the model region and, therefore, the largest 
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blockage effect. No model-induced pressure variation 
is present on the upper and lower open-jet bound- 
aries, as expected. The C p values with the slotted- 
wall calculation fall between those of the open-jet and 
solid-wall calculations, which illustrates the partially 
open behavior of the slotted wall. The slot ted- wall 
and free-air results are similar, illustrating the re- 
duction of wall interference with slotted walls. Each 
test section calculation shows a larger disturbance 
than the free- air calculation on the solid sidewall, 
although the trends among the three types of wall 
boundaries are still as expected. The calculated force 
coefficients for these cases are listed in table 2. Again, 
the trends are as expected: solid walls cause an in- 
crease in the C^, open boundaries result in a decrease 
in C^, and, slotted walls give the closest agreement 
to the free- air C p . 

Examination of Experimental Data 

Typical wall pressure measurements from the 
NTF are shown to illustrate features of the data. 
Considerable scatter has been found to exist in NTF 
wall pressure data (refs. 16 and 17). Likely, the ma- 
jority of this scatter is due to imperfections in either 
the pressure orifices or the local wall surface. Because 
scatter due to such systemic sources is generally re- 
peatable, the data quality can be greatly improved 
by subtraction of tare, or tunnel-empty, pressure val- 
ues. As noted in reference 17, the use of a tare in test 
section simulations depends on the modeling philoso- 
phy because the subtraction of a tare also eliminates 
disturbances due to other test section flow phenom- 
ena. Tared wall pressures are used in the method of 
reference 3, while the method of references 5 and 16 
uses smoothed, untared pressure data because the in- 
tent is to model the complete test section flow field. 
In the current work, the use of a tare is consis- 
tent with the problem formulation. Subtraction of 
a tare removes not only the scatter in the pressure 
data but also perturbations that result from tunnel- 
empty wall boundary layers and flow angularity, non- 
parallel walls, and other test section asymmetries 
that are not modeled. 

The NTF has nine longitudinal rows of pressure 
orifices: three in the ceiling, three in the floor, and 
three in one of the sidewalls (fig. 1). A row of 
orifices is located along the centerline of each of these 
walls. Both the ceiling and the floor contain two 
additional rows at spanwise locations of one third and 
two thirds of the test section half-width (toward the 
instrumented sidewall). Each ceiling and floor row is 
located in the middle of a solid slat that separates two 
slots. The two additional rows in the instrumented 
sidewall are located halfway between the floor and 


centerline and halfway between the centerline and 
ceiling. Typical tared and untared wall pressure 
measurements are shown in figure 6. The data were 
obtained with the Pathfinder I model at M = 0.80, 
R = 4.0 x 10 6 per foot, and Ci = 0.514. The data 
used as a tare were obtained from a tunnel calibration 
test in which the test section contained only a slender 
centerline probe that extended the length of the test 
section. The tare data were obtained at conditions 
of M — 0.80 and R — 6.1 x 10 6 per foot, the closest 
available to the Pathfinder I conditions. The data 
from the ceiling orifice rows are shown in figure 6(a). 
Data from each floor and sidewall row are shown in 
figures 6(b) and 6(c), respectively. The data show 
that the application of a tare greatly reduces the 
data scatter in all the pressure signatures. Use of the 
tare also allows identification of plugged or leaking 
orifices, and the data from such orifices are removed 
from the tared distributions. 

Note that the wall C p values shown in figure 6 
do not diminish to 0 at the test section entrance for 
either the tared or untared data. This behavior is 
a reflection of the reference values originally used 
to formulate the coefficients. Figure 6 also shows 
that the longitudinal pressure gradient on each wall 
is changed in taring the data, which indicates the ex- 
istence of tunnel-empty gradients that are subtracted 
with the tare. These gradients presumably arise from 
tunnel phenomena other than wall interference, as 
discussed above, and cannot be predicted by the cur- 
rent method. A final comment concerns the C p values 
near the outflow plane. The centerline probe begins 
to increase in diameter at station 18 to fit the NTF 
model support system. This change in diameter in- 
duces a flow gradient over the aft portion of the test 
section, which means that the wall pressures in this 
region are not true tunnel-empty values. Thus, in 
subtracting the centerline probe wall pressure data 
as a tare, the tared values downstream of station 18 
are likely to be incorrect. 

Characteristics of Test Section Solution 

Calculated and measured wall C p values arc com- 
pared in the following figures. Recall from the dis- 
cussion preceding equation (22) that the calculated 
value of C p at the test section entrance should be 
28 (because 6 is defined as C p /2). In each figure, 
the measured wall C p values have been uniformly 
shifted by a constant value so that the calculated 
and experimental C p values are in agreement at the 
test section entrance. The same offset is applied to 
each of the nine rows of pressure measurements, and 
the value of the constant depends on the value of 
8. This approach is analogous to recalculating the 
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measured wall pressure coefficients using the same 
reference conditions assumed in the simulation. 

The general behavior of the test section solution 
and the effect of varying certain simulation param- 
eters is now examined. For these illustrations, the 
geometry of the Pathfinder I model is used (fig. 4(a)). 
The model is located between tunnel stations 10.7 
and 14.9. Because the code is capable of treating only 
a single lifting surface, the modeled configuration 
consists of the wing and fuselage with no horizontal 
stabilizer. 

Effect of slot model parameters . The test sec- 
tion solution will be illustrated with the Pathfinder I 
model at a free-stream Mach number of 0.8 and an 
angle of attack of 2.70°. Other solution param- 
eters include a slot depth of 0.75 in., a slot strip 
width of 5 in., and an outflow-plane (j) xx of 0. The 
sting representation is not included. Calculated wall 
pressure distributions with two different values of 
plenum pressure are compared with measured coeffi- 
cients along the centerline of each wall in figure 7. As 
seen in figure 7(a), large discrepancies exist between 
the calculated and measured values with 5 = 0. In 
particular, the gradients upstream of the model are 
not predicted, and the model-induced acceleration on 
the ceiling is under predicted. A negative value of 5 
greatly improves the agreement in these two regions 
(fig. 7(b)), although some disagreement still exists 
between the calculated and measured C p values for 
the floor and ceiling in the region of station 10. The 
sidewall signature also shows that the location of the 
recompression aft of the wing is incorrect. In these 
figures, the decrease in pressure above the wing and 
the increase in pressure below the wing are indica- 
tive of the lift generation. Note that on the floor and 
ceiling, C p at the inflow plane is exactly equal to the 
plenum C p (that is, 25), which must be true because 
the slots begin at the inflow plane. In each case, 
the Mach number at the inflow plane is 0.80, which 
means that M 0 0 =0.80 with 5 = 0 and A/qc = 0.7964 
with 5 = -0.004. These values are determined us- 
ing equation (22). The slot volume flux q through 
each slot for these two values of 5 is shown in fig- 
ure 8. Positive q values represent outflow from the 
test section to the plenum, and negative values rep- 
resent inflow to the test section. In general, outflow 
exists along the entire length of the floor slots, and 
the magnitude of the flux decreases as the distance 
from the tunnel centerline increases. The ceiling slots 
show similar outflow up to the model location, fol- 
lowed by strong inflow through the rest of the. test 
section. Upstream of the model, the flow expands 
outward through all the slots, while the differences 
downstream of the model are due to generation of 


lift. The slots reflect the downwash aft of the wing 
by showing a general downturn of the flow in this 
region. Comparison of figure 8(a) with 8(b) shows 
that lowering the plenum pressure increases the slot 
outflow and decreases the inflow, as expected. An al- 
ternate way of viewing the effect of plenum pressure 
is to note that varying 5 changes the overall slot flux 
gradient along the slots; decreasing plenum pressure 
tends to increase the slope of q with respect to x. 
Varying the plenum pressure also affects the calcu- 
lated model lift. As 5 is changed from 0 to —0.004, 
Cl decreases from 0.5655 to 0.5601. This behavior 
makes sense qualitatively because a decrease in the 
plenum pressure aids the relieving action of the slots 
by increasing slot outflow. 

Because the cross-sectional shape of the slot edge 
is not rectangular (fig. 9), the proper slot depth spec- 
ification in the slot model is unclear. The slot lip ini- 
tially has a depth of 0.354 in., and the beveled slot 
wall continues to a depth of 1.75 in. Downstream of 
station 5 the lip is rounded, and there is a smooth 
transition to the beveled wall. A comparison of 
the wall centerline C p values obtained from solutions 
with various slot depths and the Pathfinder I model 
at M = 0.80 and a = 2.70° is shown in figure 10. Slot 
depths of 0.354 and 1.75 in. arc used as well as an 
intermediate depth of 0.75 in. The other solution pa- 
rameters are the same as in the previous case, with 
5 — -0.004. An increase in the slot depth from 0.354 
to 1.75 in. increases the model Cl by 0.002. Fig- 
ure 10 shows that increasing the slot depth increases 
the flow deceleration upstream of the model and in- 
creases the acceleration above the model. Deeper 
slots are able to maintain a larger pressure gradi- 
ent between the test section and plenum so that the 
test section appears more closed (that is, more like 
a solid- walled test section). The effective slot depth 
should be between 0.354 and 1.75 in. because the flow 
entering the slot will probably not remain attached 
along the entire depth of the beveled wall. The ef- 
fective slot depth may vary locally, depending on the 
existence (and magnitudes) of slot inflow or outflow. 
Given the shape of the slot edges, a shallow effec- 
tive slot depth is currently felt to be realistic. A slot 
depth of 0.75 in. is used in the current study. 

In the slot flow model, the effect of each slot is 
spread laterally over a slot strip of width b. The in- 
tent is to ensure a sufficient number of spanwise grid 
planes within each slot strip to adequately resolve 
the slot flow. The effect of varying b is shown in 
figures 11 and 12. The solutions are obtained with 
the Pathfinder I model at M = 0.80 and a = 2.70°, 
and the other parameters include 5 = —0.004 and 
l — 0.75 in. Figure 11(a) shows the effect of using 
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6 values from 3 to 7 in. on the wall centerline pres- 
sure coefficients. With the current grid, 6 = 7 results 
in five spanwise grid planes within each slot strip and 
seven grid planes on the solid slats between two slot 
strips. For 6 = 5 or 3, three spanwise grid planes 
fall within each slot strip. Increasing the strip width 
generally yields more positive values of C p , and in- 
creasing 6 from 3 to 5 in. causes a larger change than 
increasing 6 from 5 to 7 in. To gain further insight 
into the behavior of the slot model, figure 11(b) shows 
a comparison of the pressure coefficients that result 
when 6 = 5 and 5.5 in. Each slot strip contains five 
spanwise grid planes with 6 = 5.5 and three spanwise 
planes with 6 = 5. As shown in the figure, this change 
in strip width causes no change in the wall pressure 
coefficients. These results indicate that the solution 
is influenced more by the specified width of the slot 
strips than by the relative number of grid planes 
within the slot strips. The slot fluxes with 6 = 3 
and 7 are shown in figure 12. Increasing the strip 
width causes q to become more positive (or increases 
the gradient of the outflow). This effect is similar to 
that caused by lowering the plenum pressure. Since 
negative values of 6 are used in the current solutions, 
increasing the strip width causes the slotted wall to 
behave in a more open fashion, thus making the flux 
more sensitive to the influence of plenum pressure. 
As the strip width is increased from 3 to 5 to 7 in. at 
this angle of attack, the calculated model Cl changes 
from 0.5579 to 0.5601 to 0.5609. Because an increase 
of 6 has a diminishing effect, as shown in figure 11(a), 
a slot strip width of 7 in. is typically used in this 
study. 

Effect of test section parameters . The effects 
of the model support sting and of changing the value 
of d xx in the outflow plane are illustrated next. 
The change in the solution induced by including 
the sting representation is shown in figure 13. The 
calculated wall centerline pressures are showm for the 
Pathfinder I model at M = 0.80 and a = 2.70°. The 
sting has a greater influence on the sidewall and floor 
distributions than on the ceiling distribution for this 
lifting case. In a lifting configuration, inflow 7 exists 
from the plenum through the ceiling slots in the 
region aft of the wing, wdiich may account for the 
small effect of the sting on the ceiling. The presence 
of the sting diminishes the flow 7 velocity over the aft 
portion of the model region, as indicated by the more 
positive values of C p on the sidewall and floor. This 
increase in C p improves the agreement with measured 
values on both the sidewall and the floor. As the 
sting diameter increases downstream of the model, 
a corresponding acceleration of the flow is predicted 
on each wall. This negative C v gradient also appears 


to improve the agreement with experiment, although 
the measured values near station 18 and beyond may 
be erroneous due to the tare, as discussed above. The 
calculated Cl is 0.5601 without the sting and 0.5608 
with the sting included, wdiich shows that the sting 
has little influence on model lift. 

The effect of varying the outflow-plane value 
of (f> xx on the wall centerline C p values is showm in 
figure 14. The solution parameters are the same as 
in the previous case with the sting included. Varia- 
tion of (fr xx has no influence upstream of tunnel sta- 
tion 15. For this lifting case, fi xx has the largest effect 
on the solid sidewall and the least effect on the ceil- 
ing pressures. The small positive value of 4 > xx , which 
corresponds to a uniform acceleration of the flow 7 at 
the outflow plane, tends to improve the agreement 
with experiment (given the expected exceptions near 
the outflow plane). This variation in <j) TX causes no 
change in the calculated model force and moment 
coefficients. 

Effect of slot viscous modeling parameters . 
The effects of the viscous modeling parameters in the 
slot flow model, as implemented in reference 12, are 
now illustrated. The viscous parameters g u and rj a 
(eqs. (9) and (10)) had values of 1 in the previ- 
ous solutions presented here, which corresponds to 
a condition of no viscous losses through the slots. 
The result of independently varying the two param- 
eters is showm by considering the wall centerline pres- 
sures in figure 15. The solutions are again calcu- 
lated for the Pathfinder I model at M = 0.80 and 
a = 2.70°, with S = —0.004 and the sting represen- 
tation included. Following reference 12, 30-percent 
reductions in the viscous parameters are used. Re- 
ducing the slot width has a larger effect than reducing 
the longitudinal velocity at the slot entrance. Reduc- 
ing r] a increases the flow 7 deceleration ahead of the 
model, as indicated on each of the walls. The in- 
crease in C p continues over the entire length of the 
floor, and increased values also result aft of the wing 
(tunnel station 13) on the sidewall. The conclusion 
is that reducing the slot width makes the test section 
appear more closed, primarily in regions of slot out- 
flow. Conversely, reducing rj u results in a decreased 
flow acceleration upstream of the model and a de- 
crease in the suction peak above the wing. Rcduc- 
ing rju has almost no effect on the flow 7 at the side- 
wall. On the ceiling, the result is a leveling of the 
pressure distribution (bringing the peak values closer 
to the reference value), as expected from considera- 
tion of the small disturbance expression for pressure 
coefficient. 

The effect of a simultaneous 30-percent reduction 
in the tw 7 o viscous parameters is showm in figure 16. 
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All other solution parameters are the same as in the 
previous case. Figure 16(a) shows a comparison of 
the calculated and measured wall C p values for all 
three rows of ceiling pressure orifices. In general, 
the simulation of slot viscous effects improves the 
comparison with experiment by increasing the flow 
deceleration upstream of the model and by moving 
the location of the model-induced flow acceleration 
slightly farther downstream. However, including 
the viscous effects also decreases the magnitude of 
the suction peak in the wing region, which worsens 
the comparison with experiment. In figure 16(b), 
the results for the three rows of floor orifices are 
shown. The overall increase in Cp due to the slot 
viscous parameters results in worse agreement with 
the measured pressures. Without viscous effects, 
very good agreement exists except in the vicinity of 
stations 9-11, which is the region just ahead of the 
model. The results for the three rows of sidewall 
orifices are shown in figure 16(c). The increase in Cp 
due to viscous losses results in good agreement with 
the measured pressures along each of the orifice rows. 
The largest discrepancies are in the aft portion of the 
sting region, where the tared pressures are most likely 
to have incorrect values. 

Modifications to Slot Model 

On examination of the wall pressure coefficients 
presented so far, discrepancies between calculated 
and measured values are consistently observed in 
two regions. One area is on the ceiling aft of the 
wing-induced suction peak (tunnel station 13), and 
the other region is immediately upstream of the 
model location on both the floor and ceiling (near 
tunnel station 10). As a result of these observations, 
attempts were made to improve the simulation in 
these locations. 

The ceiling region downstream of a lifting wing is 
the only region where strong inflow from the plenum 
through the slots occurs. Because the actual plenum 
is essentially quiescent, this inflow from the plenum 
has a low momentum in comparison with the flow 
that enters the slots from the test section. In light 
of these considerations, the slot crossflow momen- 
tum equation was modified. When there is inflow 
from the plenum the quadratic term that appears in 
the momentum equation is neglected. The effect of 
excluding the quadratic term under slot inflow condi- 
tions is showm by examining the ceiling Cp values in 
figure 17. Only a slight improvement is seen in the in- 
flow region, and the effect diminishes with increasing 
distance from the centerline. Several additional ten- 
tative explanations exist for the discrepancies seen 
in this region. One possibility is that the slot flow 


model does not accurately simulate the physics of 
slot inflow regions due to either the entrance of a sig- 
nificant mass of low-energy plenum air into the test 
section or additional viscous effects within the slots. 
Another possibility is the behavior of the wind tunnel 
wall boundary layer in this region, where the adverse 
pressure gradient may cause the local wall boundary 
layer to thicken or even separate. Other causes may 
be the inability of the code to treat the blockage re- 
sulting from the model wake or flow separation from 
the aircraft and sting at angle of attack. 

Slot viscous effects may also contribute to the dis- 
crepancies seen in the region just upstream of the 
model location. Aji inflection that is not present 
in the measured pressures appears in the calculated 
floor and ceiling pressures between tunnel stations 9 
and 11; for example, refer to figures 16(a) and 16(b). 
Examination of figure 10 shows that the effect be- 
comes more pronounced as the slot depth is in- 
creased. As shown in table 1, the slots narrow to a 
minimum width at station 8.8, then gradually reopen 
to reach a constant width by station 10.6. Although 
the viscous factor % allows for a constant reduc- 
tion in slot width as implemented, viscous effects are 
likely to be more complicated in this slot-width in- 
flection region. The boundary layer on the slot wnlls 
may thicken or separate as the slot reopens, or a vor- 
tex may develop in the slot as it reopens. One result 
of these postulated models would be an effectively 
delayed reopening of the slot. As a simple means of 
examining such effects, two modified slot shapes were 
used in the test section simulation. The half- widths 
of the NTF slots as designed and of the two modi- 
fied slots are shown in figure 18 (with greatly mag- 
nified half- widths). In the first modification the re- 
opening is delayed, then occurs more gradually. The 
second modification is similar, but the minimum 
width is narrowed further. The effect of these slot 
shape modifications on the wall centerline C p values 
is illustrated in figure 19. For these calculations, r] a 
and Tj u both have values of 1. Modification 1 succeeds 
in removing the inflection between stations 9 and 1 1 
on the floor and ceiling, and good agreement with the 
measured pressures exists at all other locations. Ad- 
ditionally, modification 2 improves prediction of the 
deceleration near station 10. Although the results of 
these modifications are encouraging, they offer little 
insight into the actual flow' physics in such regions. 
The effects of these slot modifications are examined 
further in the next section. 

Simulation Performance 

In this section, the parameters of the simula- 
tion are adjusted to give the best possible agreement 
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with the experiment. The Pathfinder I and Boe- 
ing 767 models are considered at both zero lift and 
cruise lift, and the comparisons are carried out at 
a Mach number of 0.80 to approximate the cruise 
conditions of the aircraft. At the higher lift coeffi- 
cients {Ci > 0.5), significant regions of supersonic 
flow exist at this Mach number. The experimen- 
tal data with the Pathfinder I model were obtained 
at R = 4.0 x 10 6 per foot, and the Boeing 767 data 
were obtained at R = 7.5 x 10 6 per foot. Although 
almost a factor-of-2 difference exists, both values are 
low with respect to full-scale flight values and to the 
NTF operating envelope. The same tare data, ob- 
tained at M = 0.80 and R — 6.1 x 10 6 per foot, are 
applied to each set of measured wall pressure data. 

Pathfinder I. The case that has been used for 
illustration so far, the Pathfinder I model at cruise 
lift, is considered first. These calculations are per- 
formed with M = 0.80, outflow cf> xx = 0.0001, and 
the sting representation included. The slot param- 
eters are 6 — —0.004, Slot depth = 0.75 in., and 
Slot strip width = 7 in. Slot shape modification 1 
is used unless otherwise noted. With an angle of 
attack of 2.70°, as used in the previous cases, the 
calculated model Ci is 0.562. To match the experi- 
mental value of 0.514, an angle of attack of 2.30° is 
used here. The wall centerline pressure coefficients 
resulting from these two angles of attack are com- 
pared in figure 20. Slight differences are seen through 
the model region, and the most significant effect is a 
change in the C p level on the floor. An angle of attack 
of 2.30° is used for all other calculations shown for 
this case. Figure 21 shows the result of including slot 
viscous losses with 15-percent reductions for both rj a 
and rj u . Very good agreement with experiment is 
seen. The inclusion of viscous effects improves the 
prediction of the suction peak in the model region on 
both the sidewall and ceiling. However, the pressures 
on the floor are in better agreement with no viscous 
losses. This observation raises the possibility that 
the effects of viscosity on the slot flow vary, depend- 
ing on the local condition of inflow or outflow — above 
or below a lifting wing, for example. Including the 
viscous losses causes the calculated Ci to increase 
from 0.5145 to 0.5152, Cj \j to change from —0.1515 
to —0.1511, and C]j to increase by 0.0002. These 
increments are quite small. The solution illustrated 
in figure 21, where rj a = rj u = 0.85, is considered the 
best overall prediction for the current case. 

A more detailed examination of the effects of the 
slot width modifications is informative. In figure 22, 
the wall pressures from all nine rows of pressure ori- 
fices are shown. Results are shown for the actual slot 
geometry as well as for both modified slot shapes, and 


each calculation is performed with rj a = T) u = 0.85. 
The case with slot modification 1 is the same as that 
presented in figure 21. Examination of the ceiling 
pressures (fig. 22(a)) shows that on the wall center- 
line the modified slot shapes give better agreement 
with experiment between stations 7 and 11, and mod- 
ification 2 results in the best agreement. On the other 
two rows, the actual slot shape and modification 1 
tend to give better agreement in this region. Al- 
though part of the uncertainty is due to the scarcity 
of data in the region, the values on the centerline 
may differ for another reason. In this region, the 
flow is decelerating due to the blockage imposed by 
the aircraft model, and the blockage induced by the 
fuselage should influence the pressure along the wall 
centerline more than at the outer locations. Since 
the fuselage begins at station 10.7, it is probable that 
the blockage induced by the fuselage is not properly 
captured by the small disturbance numerical formu- 
lation. The three floor pressure rows are shown in 
figure 22(b), and the solutions with slot modifica- 
tions 1 and 2 are both adequate. Along the three 
sidewall pressure rows (fig. 22(c)), good agreement 
with experiment is seen when modification 1 is used. 
Overall, consideration of data from all the orifice 
rows instead of just the wall centerline rows is clearly 
helpful when making comparisons between measured 
and computed values. This approach allows span- 
wise and vertical variations to be considered and al- 
lows scarcity of data in some regions to be overcome. 
Finally, note that the slot modifications have no ef- 
fect on the previously discussed discrepancies on the 
ceiling aft of the suction peak (fig. 22(a)). 

The Pathfinder I model at zero lift is considered 
next. All the slot model parameters have the same 
values that resulted in the best agreement with ex- 
periment in the cruise-lift case. The calculated wall 
Cp values at an angle of attack of —2.0° are compared 
with measured values in figure 23. Calculations are 
shown with both the design slot shape and slot mod- 
ification 1. The slot width modification has a pro- 
nounced effect on the calculated pressures through- 
out the model region, which results in an improved 
agreement with measured values that is similar to 
the agreement shown in the lifting case. The cal- 
culated lift coefficients are 0.0001 with the design 
slot and 0.0021 with the modified slot, while the Ci 
corresponding to the measured pressures is 0.011. 
Comparison of figure 23(a) with 23(b) shows that 
the calculated ceiling and floor distributions are al- 
most identical, as expected from a zero-lift calcula- 
tion. The measured floor pressures, however, exhibit 
a slightly higher suction through the model region 
than the ceiling pressures exhibit. Such behavior 
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would be expected if the model were generating a 
slight negative lift, but the measured Cl is positive. 
This discrepancy is a result of the taring of the pres- 
sure data because empty-tunnel gradients are sub- 
tracted as part of the tare. Also, the asymmetry 
in the measured pressures may be partially due to 
the influence of the horizontal stabilizer, which is not 
modeled in the calculation. Excellent agreement be- 
tween measured and calculated values exists along 
the ceiling pressure rows with the modified slot shape 
(fig. 23(a)), but on the floor the low pressure through 
the model region is underpredicted (fig. 23(b)). Fig- 
ure 23(c) shows excellent agreement in the sidewall 
pressures with the modified slot shape, with the ex- 
pected exception of those pressures near the outflow 
plane. The proper prediction of the C p on the solid 
sidewall is an indication that the openness ratio of the 
slotted walls that exists experimentally is properly 
simulated in the calculation. 

Boeing 767, The Boeing 767 model at zero lift 
is considered next. The 767 model is slightly larger 
than the Pathfinder I model and has a fuselage ex- 
tending from tunnel station 10.8 to station 15.4. The 
values of the slot model parameters used here are 
the same as those that give the best predictions with 
the Pathfinder I model. A higher value of 4>xx in 
the outflow plane is necessary for the Boeing 767 
model due to the different sting arrangement. A 
value of 0.0003 was found to give the best behav- 
ior and is used for the results showm here. Figure 24 
shows the wall pressures that result from calculations 
performed both with and without slot-width modifi- 
cation 1. Calculated Cl values are 0.0045 with the 
design slot and 0.0071 with the modified slot, while 
the measured Cl is 0.028 for this case. On the ceil- 
ing, figure 24(a) shows that the solution with the un- 
modified slot shape gives good agreement with the 
measured pressures. Again, the largest discrepancy 
is on the centerline at stations 8 and 9. Both solu- 
tions tend to overpredict the negative pressures ex- 
isting through the model region. Ahead of the model, 
the comparison on the floor is similar to that on the 
ceiling (fig. 24(b)). Through the model region the 
solutions with the two different slot shapes bracket 
the measured values, while in the sting region the C p 
values are undcrpredicted though the trends are cor- 
rect. Comparison of figure 24(a) with 24(b) shows 
that the calculated pressures again are similar on 
the floor and ceiling, while the measured pressures 
show a slight asymmetry through the model region 
and a significant asymmetry through the sting region. 
In the sting region, these observations are consistent 
because the Boeing 767 sting is located above the 
fuselage (fig. 4(b)), whereas in the code the fuselage 


and sting are located on the same axis. The side- 
wall pressures (fig. 24(c)) support the observations 
made on the ceiling and floor, which include the over- 
prediction of the pressures in the model region in the 
upper half of the test section, the underprediction 
in the lower half, and the asymmetry in the sting 
region. As in the case of the Pathfinder I model 
near zero lift, the measured pressures in the model 
region have a characteristic that would result from 
the model generating a slight negative lift: higher 
acceleration is found on the floor than on the ceiling. 
Because the measured lift is actually a positive value, 
the discrepancy between measured floor and ceiling 
pressures is again attributed to the tare subtraction. 
The current solution with the design slot shape gives 
good overall agreement with the measured pressures, 
but these comments indicate that improved agree- 
ment with the measured pressures might be achieved 
by decreasing the angle of attack in the calculation. 
Such a change is analogous to correction for tunnel- 
empty flow angularity. Figure 25 shows the effect 
of lowering the angle of attack by 0.3°. The calcu- 
lated Cl is —0.0328 at this angle of attack. Improved 
agreement with the measured pressures through the 
model region appears along all the orifice rows. 

Finally, consider the Boeing 767 model at cruise 
conditions. The wall centerline C p values resulting 
from a solution at an angle of attack of 2.15° with 
the design slot shape and the same values of the slot 
model parameters used in the previous cases are com- 
pared with measured values in figure 26. The calcu- 
lation predicts an incorrect location for the begin- 
ning of the suction peak above the model and for 
the recompression aft of the model on the sidewall. 
The C p level on the floor in the model region is un- 
derpredicted. These types of behavior are improved 
by lowering the plenum pressure coefficient. The ef- 
fect of decreasing 6 to —0.005 while all other slot 
parameters remain unchanged is shown in figure 27. 
Figure 27 also shows the effect of including slot mod- 
ification 1. The calculated Cl — 0.5464 with the de- 
sign slot and 0.5479 with slot modification 1, and the 
measured Cl = 0.548. On the ceiling (fig. 27(a)), 
the solution with the unmodified slot shape gives 
good agreement with measured values up to the suc- 
tion peak over the wing. The disagreement down- 
stream of the suction peak is largest on the center- 
line and diminishes toward the sidewall. Given 
this spanwisc variation, much of the discrepancy 
is attributable to improper simulation of the high- 
mounted sting (including sting/fuselage interference 
and wake/separation effects). The behavior near 
tunnel station 10 is similar to that observed in the 
previous cases. The design slot shape also gives 
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good results on the floor pressure rows, as shown 
in figure 27(b). Again, the largest disagreement is 
found near station 10. Very good agreement be- 
tween measured and calulated C p values exists on 
the sidewall when the unmodified slot shape is used 
(fig. 27(c)). 

Discussion of results . The goal of the current 
work is to derive a treatment of the slotted-wall flow 
field that allows consistent prediction of the test 
section flow. Because of the many parameters in 
the current slotted-wall model, a proper combination 
is difficult to select. A combination of parameter 
values is developed here that gives reasonably good 
predictions for two models at lifting and nonlifting 
conditions (although in the final case the simulation 
is improved by further lowering the plenum pressure). 
These results are encouraging, but some additional 
issues are identified. 

The results of the slot width modification shown 
here indicate that a delayed reopening of the slots 
downstream of their minimum in width, such as 
might result from the effects of viscosity, improves 
the agreement with measured wall Cp values. With 
the Pathfinder I data, the slot modification tends to 
properly compensate the calculation in this region. 
With the Boeing 767 data, the discrepancy is smaller, 
and the same slot modification appears to over- 
compensate in this region. Although no mechanisms 
for the width reduction have been identified, the 
simple approach taken here shows that the effect 
is independent of lift. The behavior is tentatively 
attributed to the test section wall boundary layers 
and/or to viscous effects in the slots. Because the 
Reynolds number is lower in the Pathfinder I data 
than in the Boeing 767 data, larger viscous effects 
in the slots can be expected with the Pathfinder I 
data. The effects of viscosity on the walls and in 
the slots are further clouded by the subtraction of 
tunnel-empty wall C p values as a tare. 

The subtraction of tares from the measured wall 
pressures exposes concerns that pertain to wall inter- 
ference corrections. As evident from the zero-lift 
comparisons, an offset between the measured and cal- 
culated Cl values results when the calculated wall 
C p values are brought into agreement with the tared, 
measured values. This tare- induced offset is indica- 
tive of global flow angularity and Mach number off- 
sets that are discarded with the tare. Such offsets 
would have to be considered, perhaps through anal- 
ysis of the actual tunnel-empty wall pressure data, 
before attempting to use the code results to calculate 
any actual corrections to wind tunnel data. 


The outflow-plane <j> xx specification in the current 
method is intended to mimic the gross C p trends 
near the reentry region. It is apparent from the 
wall C p results that this approach gives reasonable 
results with the Pathfinder I sting arrangement but 
less accurate results with the high-mounted sting of 
the Boeing 767. Although the calculated model coef- 
ficients seem insensitive to the outflow specification, 
they may be affected by proper sting representation 
because the sting can appreciably influence the flow 
gradients in the aft part of the model region. 

Finally, the current results must be viewed in the 
proper context. The simulation method is capable of 
reproducing the wall pressure signature that exists 
in the NTF test section. An implicit assumption is 
that the effect of the walls on the model flow field is 
similarly predicted. Although a relationship between 
the wall effects and the model flow field certainly 
exists, this assumption remains unproven due to the 
lack of appropriate experimental data from the model 
surface. In the following section, the calculated effect 
of the test section on the model flow field is shown. 

Wall Interference Illustration 

The effect of the test section on the aircraft flow 
field is shown by comparing the results of free-air and 
test section calculations with the Pathfinder I model 
at M = 0.80 and an angle of attack of 2.70°. Test 
section simulations are performed with the sting in- 
cluded and slot shape modification 1 used for both 
discrete and homogeneous slotted-wall boundaries. 
The discrete slotted-wall parameters have the opti- 
mum values that were developed in the previous sec- 
tion. Table 3 lists the calculated force and moment 
coefficients and shows that the homogeneous calcu- 
lation gives results that are almost identical to the 
free-air values. The discrete slotted-wall calculation 
shows an underprediction of the lift by 0.7 percent, 
of the pitching moment by 1.5 percent, and of the 
drag by 7.7 percent in comparison with the free-air 
calculation. 

Local variations induced by the test section are 
examined by considering the wing pressure distribu- 
tion. For the free-air and discrete slotted-wall calcu- 
lations, pressure distributions at wing semispan lo- 
cations of 21 percent, 52 percent, and 82 percent 
are shown in figure 28. These comparisons show 
that the shock locations predicted by the two cal- 
culations differ by less than 1 percent of the chord at 
each spanwise location. The test section calculation 
shows a greater upper surface expansion at the in- 
board and midspan locations (figs. 28(a) and 28(b)), 
which indicates the existence of slightly stronger 
shock waves. The outboard station shows the only 
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substantial difference between the lower surface pres- 
sures (fig. 28(c)), which probably illustrates a flow 
angularity effect. Although small, the differences 
shown here are likely to increase with models larger 
than the Pathfinder I. The calculated wing pres- 
sure coefficients with frec-air and homogeneous slot- 
ted boundaries are identical; the 52-percent semispan 
location shown in figure 29 is an example. 

As shown in the preceding sections, good agree- 
ment with measured wall pressures is obtained by 
using negative plenum pressure coefficients in the 
discrete-slot model. Reference 15 shows that the 
discrete-slot model with <5 = 0 gives results very sim- 
ilar to those obtained using the homogeneous slotted 
boundary. (Also, compare figs. 5 and 7(a).) Thus, 
a more accurate portrayal of wall interference effects 
is probably given by the current discrete slotted-wall 
calculation. 

Concluding Remarks 

A computational method was developed that is 
capable of simulating the flow in a transonic test 
section with longitudinally slotted walls. The slotted 
walls may be modeled either with discrete slots or in 
the classical homogeneous manner. The simulation 
results were compared with measured wall pressure 
data to judge the performance of the method at the 
walls. For the aircraft model geometries and limited 
available data shown here, the discrete-slot method 
was found to give consistently good results with 
minimal variation of the parameters that appeared in 
the slot flow model. The value of the plenum pressure 
coefficient was one of these parameters, and negative 
values (on the order of -0.01) were required to yield 
proper behavior. Because the classical homogeneous 
formulation does not have a means to vary plenum 
pressure, this technique was incapable of properly 
predicting the measured wall pressures. 

The simulation method may be used to predict 
the influence of the test section environment on the 
aircraft flow field. These predictions are made by per- 
forming both free-air and test section calculations, 


then examining the incremental differences between 
the two solutions. Such increments may then be 
considered in conjunction with data measured dur- 
ing a wind tunnel test. These increments cannot be 
used as wall interference corrections at present due to 
the flow angularity uncertainties that arise from the 
subtraction of tunnel-empty wall pressures as a tare. 
However, the simulation can give indications of the 
variations expected in force and moment coefficients 
as well as in measured model pressure distributions. 

Additional wall pressure data are needed to eval- 
uate the simulation method more thoroughly. Data 
obtained with simple wing/ fuselage model configura- 
tions at high Reynolds numbers would be particularly 
useful. The Reynolds numbers of the data sets used 
here are on the low end of the NTF operating en- 
velope, and how the simulation will perform or how 
the slot model parameters might need to be changed 
when confronted with high Reynolds number data 
(up to an order of magnitude larger than considered 
here) is not known. Some localized disagreement be- 
tween the measured and calculated wall pressure co- 
efficients is currently noted that may result from im- 
proper modeling of slot and/or wall viscous effects. 
Because these viscous issues are complicated by the 
subtraction of tunnel-empty tare pressures, a clearer 
determination requires examination of data at higher 
Reynolds numbers. 

The ultimate intent is to use the current method 
as a pretest indicator, showing whether wall interfer- 
ence is a concern for a given model and flow condi- 
tion. The ability of the method to simulate the flow 
at the test section walls has been shown here. Due 
to the model sizes, very little wall interference exists 
with the models considered here. Thus, the useful- 
ness of the method to predict wall interference effects 
on the model requires validation with data that do 
not currently exist. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
June 10, 1993 
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Appendix A 

Solution of Centerplane Slot Flow 
Equations 

The integration of the two slot flow equations is 
not a straightforward procedure due to the complex- 
ity of the expressions. In the auxiliary equations, the 
variables of interest cannot be isolated in a closed 
form. For convenience, the centerplane slot equations 
(eqs. (4), (6), (18), (19), and (20)) are summarized 
and repeated as follows: 
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The goal is to simultaneously integrate equations (Al) 
and (A2) along the slot, given an estimate of q s . 
However, the auxiliary relations given by equa- 
tions (A3) and (A4) are necessary to define the 
terms Q and V. These relations are expressed in 
terms of the transform variable s, while the differen- 
tial equations are in terms of the variable y p . Equa- 
tion (A 5) gives the relationship between y p and s, 
but s cannot be isolated in this expression. An it- 
erative process is necessary to calculate s, given y p . 
An alternate approach is to rewrite equations (Al) 
and (A2) in terms of s. This approach is simple 



for equation (A2) and requires only direct substitu- 
tion of equations (A3), (A4), and (A5). To rewrite 
equation (Al), however, dy p /dx must be replaced 
by ds/dx. A derivation of this relation is now pre- 
sented. 

Noting that equation (A5) expresses y p /a s as a 
function of 5, it is convenient to write the simple 
expression 
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Using equation (A5) it can be shown that 
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and substitution of equation (A7) into equation (A6) 
yields the relation 
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Equation (A8) is the desired relation between dy p /dx 
and ds/dx. Substituting equation (A8) into equa- 
tion (Al) and rearranging yields 
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All terms in equation (A9) are either known functions 
of the independent variable x or explicit functions 
of the dependent variable s . Finally, recall that 
equation (5) defines the location of the limiting value 
for y P) denoted by y p>0 . Because the variable of 
integration is now s instead of y p , equation (A7) is 
used to obtain the value of s that corresponds to y PjQ . 
Let this limiting value of s be denoted by s p , a . During 
the integration of equations (A9) and (A2), s must be 
reset to the value s P:Q whenever the projected value 
of s would exceed that allowed by s Pi0 . 

To summarize, the differential equation for dy p /dx 
has been replaced by an equivalent differential equa- 
tion for ds/dx. The result is a system of two differen- 
tial equations, equations (A9) and (A2), that contain 
the two unknown functions 5 and p. A simultane- 
ous integration of the two equations along each slot 
yields s and p. 
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Appendix B 

Global Coarse Grid Generation 

This section details the generation of the global 
finite difference grid. The physical domain described 
by the Cartesian coordinates (x, y , z) is related to 
the computational domain given by The 

longitudinal direction is given by x and £, the span- 
wise direction by y and 77 , and the vertical direction 
by z and £. The grid spacing in each coordinate 
direction is constant in the computational domain, 
which greatly simplifies the finite-difference formu- 
lation. The transformations relating the two do- 
mains are simple algebraic expressions. To avoid sin- 
gularities in the transformation, the transformation 
derivatives (metrics) to second order are smooth and 
continuous. 


In the NTF test section, station 13 (13 ft from 
the beginning of the test section) is the center of 
sting rotation. The model wing is typically near this 
position. The grid in the x direction is symmetric and 
is centered around NTF station 13. For convenience, 
£ extends from -1 to 1 , and the corresponding 
x values cover negative to positive infinity. In the 
x direction, the physical domain is divided into five 
regions, with the middle region centered on station 13 
(fig. Bl). Because the grid is symmetric, the two 
zones ahead of the wing region are treated the same 
as the two zones dowmstream of the wing. The central 
region (zone 1 ) contains the wing and may be thought 
of as extending from the leading edge of the root to 
the trailing edge of the tip (for a typical swept wing). 
This region is slightly expanded as necessary to locate 
the midpoint of the region at station 13. Sixteen 
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Figure Bl. Zones used in grid generation in vertical ( z ) and streamwise (x) directions. 


cells are uniformly distributed in zone 1. The trans- 
form is written as 

x = A 0 + Ait (Bl) 

Because £ = 0 at station 13 and the values of x 
and £ at the edge of the wing region are known, the 
constants A$ and A\ arc easily evaluated. The next 
zones cover the regions that extend from the edges 
of the wing zone to the test section inflow (zone 2 
upstream) and outflow (zone 2 downstream) planes. 
(Sec fig. Bl.) In these regions, a cubic function is 
used to give a mild stretching yet still maintain a 
reasonable spacing at the test section boundaries. 
The transform is given by 

x — Aq + + A‘2 QLA) (B2) 

where £i is the value of £ at the edge of the wing 
region and £2 is equal to the value of £ at infinity 
(that is, £2 = — 1 upstream and +1 downstream). 
The constant A^ is determined by the locations of 
the inflow and outflow planes. Note that as the edges 
of zone 1 are approached from these regions, £ — > £j, 
and equation (B2) reduces to equation (Bl). Using 
superposition in this way ensures that the metrics are 
smooth and continuous at the boundaries between 
the regions. The grid-generation functions in both 
the x and y directions arc built up in this manner. 
The final longitudinal regions, defined by zone 3, 
extend from the inflow plane to upstream infinity and 
from the outflow plane to downstream infinity. The 
grid-generation function in these regions is written as 

x = ^ + ^ + 42 (&=% )’ - + Ai *“ [H^f )*] (B3) 

where £3 is the value of £ at the test section boundary. 
The constant ^3 is evaluated by requiring the next- 
to-last grid point to be located a distance Aq from the 
boundary (that is, one half the test section length). 
Note that as £ approaches its maximum value £2, the 
argument of the tangent function approaches 7 t/ 2, 
which means that the function itself tends toward 
infinity. 

In the spanwise direction, 77 ranges from a value 
of 0 at the tunnel centerline to a value of 1 at 
positive infinity. Constant spacing is used between 
the centerline and the test section wall. The spacing 
is chosen to align a grid plane with the slot closest 
to the tunnel centerline. Due to the geometry of the 
NTF test section, this technique ensures that grid 
planes are also aligned with the other two slots as well 


as with the test section sidewall. The transformation 
equation is written as 

V = BiV (B4) 

where B\ is determined by specifying the number 
of grid planes that lie between the centerline and 
the first slot. Outside the test section, the tangent 
function is again used to stretch the grid toward 
infinity. The expression is 

y = B\t] + B<2 tan ^ (— — — ) (B5) 

2 \m-mJ 

where 771 is the value of 77 at the wall, and 772 is the 
maximum value of 77. The constant B<i is evaluated 
by requiring that the next-to-last grid point is located 
a distance from the centerline that is equal to twice 
the y location of the wall. 

The grid generation in the z direction is compli- 
cated by two factors. One constraint is the desire to 
have a grid symmetric about the wing plane. How- 
ever, the wing plane may not be equidistant from the 
floor and ceiling. The other factor is the need to min- 
imize the total number of grid points while maintain- 
ing adequate grid resolution near both the wing and 
the walls. These demands arc accommodated by di- 
viding the region inside the test section into three re- 
gions (fig. Bl): a central symmetric region (zone 1), a 
region below the ceiling (zone 2 upper), and a region 
above the floor (zone 2 lower). The same transfor- 
mation is used near the floor and the ceiling, but the 
constants in the transform equation differ if the wing 
plane is not centered between the floor and ceiling. 
The test section centerline is located at z = 0, and 
the wing plane is at £ = 0. Thus, if the wing plane 
is located at z — 0, then the entire vertical grid is 
symmetric about z = 0. The ceiling and floor are at 
£ locations of ±1. Zone 1 encompasses half the test 
section height and is centered about the wing plane. 
The transformation in this region is written as 

x = z w + Cl tan (|c) (B6) 

where z w denotes the position of the wing plane. The 
constant C\ is evaluated by specifying the number of 
grid points that lie in this region. The regions defined 
by zone 2 extend from the edges of the central region 
to the floor and ceiling. The transformation is given 
by 

Z = *w + C 2 tan [| (C - Cw)] (1 + DiO (B7) 
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where the subscript W denotes the value at the wall. 
This transform serves to slightly refine the spacing 
in the vicinity of the wall. Equation (B 7 ) has the 
same form as equation (B6); both are simple tangent 
functions, but equation (B6) is centered on the wing 
plane while the tangent function in equation (B 7 ) 
is shifted so that the function is centered on the 
wall. The term containing the constant D\ is a 
scaling factor that can be thought of as expanding or 
compressing the grid in these regions. To illustrate, 
recall that equation (B 7 ) applies to both the floor 
and ceiling regions, where the constants C2 and D\ 
have the same values for these regions only if the wing 
plane is centered between the floor and ceiling. Let 
the resulting grid be thought of as uncompressed. If 
the wing plane is located at some negative z value, 
then the (1 + D\Q term will serve to compress the 


grid near the floor and expand the grid near the 
ceiling, and C2 and D\ will have different values 
in the two regions. Solving for the two constants 
requires two constraints. The constraints are that 
equations (B6) and (B 7 ) have the same value and the 
same first derivative at the edge of zone 1 . Outside 
the test section walls (zone 3 ), a tangent function 
is used to stretch the grid toward infinity. The 
transformation is written as 


z = z\y + C3 tan 


7 r 
2 


( C-Civ 

VC2 - C w 


(B8) 


where C2 is the value of C corresponding to z = 00. 
The constant C3 is 1 evaluated by requiring the first 
derivatives of equations (B 7 ) and (B8) to match at 
the wall. 
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Table 1. Ordinates Defining NTF Slot Shape 

Ordinates given as slot half-widths 
measured from slot centerline 


Tunnel 

station, 

ft 

Slot 

half-width, 

in. 

0.604 

0.000 

.686 

.010 

.754 

.032 

.912 

.100 

1.090 

.176 

1.281 

.257 

1.481 

.338 

1.689 

.413 

1.907 

.487 

2.138 

.552 

2.384 

.611 

2.651 

.663 

2.944 

.705 

3.072 

.720 

3.484 

.751 

4.144 

.757 

4.695 

.729 

5.176 

.691 

5.612 

.652 

6.015 

.612 

6.401 

.567 


Tunnel 

Slot 

station, 

half-width, 

ft 

in. 

6.776 

0.518 

7.149 

.470 

7.529 

.417 

7.923 

.367 

8.150 

.341 

8.314 

.324 

8.478 

.309 

8.642 

.298 

8.806 

.293 

8.970 

.294 

9.134 

.301 

9.298 

.316 

9.463 

.336 

9.627 

.361 

9.791 

.389 

9.954 

.417 

10.118 

.444 

10.282 

.467 

10.446 

.484 

10.611 

.492 

19.667 

.492 



Table 2. Effect of Homogeneous Outer Boundary Conditions 
on Calculated Pathfinder I Force and Moment Coefficients 

[M = 0.80 and a = 2.70°] 


Boundary 

C L 

Cm 

C D 

Free air 

0.5669 

-0.1573 

0.0183 

Open jet 

.5618 

-.1579 

.0182 

Solid 

.5772 

-.1570 

.0187 

Slotted 

.5659 

-.1578 

.0183 


Table 3. Comparison of Calculated Force and Moment Coefficients 
With Outer Boundary Conditions 

Pathfinder I at M = 0.80 and a = 2.70°; sting] 
included in test section calculations j 


Boundary 

C L 

Cm 

Co 

Free air 

0.5669 

-0.1573 

0.0183 

Homogeneous slotted-wall 

.5667 

-.1566 

.0184 

Discrete slotted-wall 

.5632 

-.1549 

.0169 
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Figure 1. NTF test section. Dashed lines indicate locations of wall pressure orifice rows. All dimensions are 
in feet. 
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Figure 2 . Slot flow model in slot centerplane. 
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Figure 3. Cross section of slot model geometry for isolated slot showing specification of boundary conditions. 
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Figure 9. Sketch of NTF slot lip shape. All linear dimensions are in inches. 
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Figure 18. Slot width modifications 1 and 2 compared with actual NTF slot shape. 
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